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Abstract 

Subsurface  radar  was  used  to  profile  ice  and  snow  conditions  on  the  Ross  Ice 
Shelf  at  McMurdo  Station,  Antarctica,  during  mid-January  1 993.  Deconvolution 
and  migration  were  often  used  to  improve  vertical  resolution  and  spatial 
imaging.  Profiles  ata  pulse  centerfrequency  of 400  MHz  along  the  3.2-km-long 
Pegasus  ice  runway  show  many  low-density  horizons  above  9  m  depth  that  are 
up  to  30  m  long.  They  are  associated  with  air  bubbles  included  during  refreezing 
of  meltwaterand  are  interpreted  as  layers  between  a  few  and  tens  of  centimeters 
thick.  There  is  a  strong  reflecting  horizon  atabout9  m  depth  that  is  probably  from 
brine  intrusion  as  it  is  continuous  with  the  intrusion  into  the  snow  to  the  east. 
Diffraction  asymptotes  give  a  dielectric  constant  near  3.2  for  material  above  the 
brine  level,  a  value  that  implies  near-solid  ice,  Profiles  at  TOO  MHz  along  the 
road  between  Pegasus  runway  and  Williams  Field  in  the  accumulation  zone 
show  snow  features  such  as  layer  deformation  and  intrusive  brine  layers  that 
both  abruptly  and  gradually  change  in  depth,  A  single  profile  at  a  relic  solid 
waste  dump  at  Williams  Field  detected  buried  debris  and  ice  within  the  upper 
7  m.  A  survey  of  a  suspected  fuel  spill  shows  some  local  disturbances  nearthe 
center,  but  no  excavation  was  done  to  verify  the  findings.  Profiles  traversing  the 
sewage  sumps  at  Williams  Field  outline  the  extent  of  the  sewage  deposition,  and 
give  depths  to  contaminated  snow  that  closely  agree  with  observation.  Despite 
variability  in  dielectric  properties,  single-layer  migration  effectively  improves  the 
resolution  of  subsurface  conditions.  Recommendations  are  made  for  future 
surveys. 


Cover:  Radar  profiling  at  400  MHz  pulse  center  frequency  along  the  Pegasus 
runway,  January  1993. 
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PREFACE 


This  report  was  prepared  by  Dr.  Steven  A.  Arcone,  Geophysicist,  Allan  J.  Delaney, 
Physical  Science  Technician,  both  of  the  Snow  and  Ice  Branch,  Research  Division,  and  W ayne 
Tobiasson,  Civil  Engineer,  of  the  Civil  and  Geotechnical  Engineering  Research  Branch, 
Experimental  Engineering  Division,  U.S.  Army  Cold  Regions  Research  and  Engineering 
Laboratory.  Funding  for  this  research  was  provided  by  the  National  Science  Foundation, 
Grant  no.  DPP-872  0063,  and  the  Office  of  the  Chief  of  Engineers  through  DA  Project 
4A161102AT24,  Research  in  Snow,  Ice  and  Frozen  Ground ;  Task  SS,  Properties  of  Cold  Regions 
Materials;  Work  Unit  E10,  Electromagnetic  Characteristics  and  Dynamics  of  Snow,  Ice  and  Frozen 
Ground. 

This  report  was  technically  reviewed  by  Dr.  Anthony  Gow  and  George  Blaisdell  of 
CRREL. 

The  contents  of  this  report  are  not  to  be  used  for  advertising  or  promotional  purposes. 
Citation  of  brand  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use 
of  such  commercial  products. 
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Subsurface  Radar  Investigations  at  the 
Pegasus  Glacial-Ice  Runway  and 
Williams  Field,  McMurdo  Station,  Antarctica 

STEVEN  A.  ARCONE,  ALLAN  J.  DELANEY  AND  WAYNE  TOBIASSON 


INTRODUCTION 

This  report  discusses  the  use  of  Ground-Pen¬ 
etrating  Radar  (GPR)  to  detect  and  map  artificial 
and  natural  subsurface  conditions  at  the  Pegasus 
Glacial-Ice  Runway  (henceforth  referred  to  as  Pe¬ 
gasus  runway),  the  access  road  to  Pegasus  runway 
and  the  Williams  Field  support  facilities  on  the  Ross 
Ice  Shelf  near  McMurdo  Station  in  Antarctica. 
GPR  has  been  a  common  geophysical  tool  for 
near-surface  geologic  exploration  in  cold  regions 
(e.g.,  Annan  and  Davis  1976,  Kovacs  and  Morey 
1979,  Arcone  and  Delaney  1987,  Delaney  et  al. 

1991,  Arcone  et  al.  1992,  Arcone  et  al.,  in  prep.  ') 
owing  to  the  excellent  penetration  of  radio  waves  in 
ice  and  frozen  soils  and  the  vertical  resolution  of  the 
short  pulse.  GPR,  operating  at  a  variety  of  band- 
widths,  can  delineate  subsurface  features  in  snow 
and  ice  at  depths  ranging  from  less  than  1  to  over 
100  m.  Vertical  resolution  is  determined  by  the 
bandwidth  center  frequency  of  the  pulse  used,  and 
is  generally  about  20  cm  in  ice  for  antennas  operat¬ 
ing  near  500  MHz,  and  less  for  materials  with  a 
higher  dielectric  constant  such  as  permafrost.  This 
exceptional  resolution  has  also  led  to  GPR  becoming 
a  popular  device  for  utility  detection  and  toxic  and 
hazardous  waste  site  characterization  (e.g.,  Pilon 

1992,  Hanninen  and  Autio  1992)  because  subsurface 
debris,  water  and  fuel  spills,  or  the  soil  disturbances 
they  cause,  can  present  a  strong  radar  target. 

The  objectives  of  these  surveys  were  to  character¬ 
ize  the  ice  and  snow  structure  above  the  brine  layer 
in  the  Pegasus-Williams  Field  area  to  support 
planned  and  ongoing  operations  projects,  and  to 
evaluate  the  ability  of  GPR  to  locate  fuel  spills,  solid 
waste  dumps  and  other  conditions  that  might  re¬ 
quire  remediation.  GPR  transducers  operating  at 
both  100-  and  400-MHz  pulse  center  frequencies 


were  towed  by  truck  along  lengthy  profiles  and  by 
hand  along  short  transects.  Signal  processing  using 
deconvolution  and  migration  was  used  to  improve 
vertical  resolution  and  spatial  imaging.  Limited  cor¬ 
ing,  probing  and  excavation  verified  interpretation 
of  GPR  features  of  interest. 

Radio  echo  sounding  and  GPR  have  been  used 
extensively  on  the  Ross  Ice  Shelf  in  the  McMurdo 
Station  vicinity  (Clough  1973;  Kovacs  and  Gow 
1975,  1977a,b;  Kovacs  et  al.  1982)  and  in  other  ice 
shelf  locations  far  removed  from  McMurdo  (e.g., 
Jezek  1980,  Morey  and  Kovacs  1982).  The  most  ex¬ 
tensive  studies  near  McMurdo  were  by  Kovacs  and 
Gow  (1975)  and  Kovacs  et  al.  (1982),  who  included 
drilling  and  ice  core  analysis  with  their  radar  pro¬ 
files.  Their  work  concentrated  on  characterizing  the 
brine  layer  (an  intrusive  layer,  probably  of  finite 
thickness,  and  fed  mainly  at  the  front  of  the  ice 
shelf),  reflections  from  which  are  the  dominant  fea¬ 
ture  in  their  records.  Their  data  contain  responses  to 
many  of  the  natural  features  such  as  cracking,  fold¬ 
ing  and  abrupt  changes  in  brine  level.  Responses  to 
these  features  are  presented  here  also,  but  with  the 
advantages  of  a  higher  resolution  bandwidth.  We 
also  use  digital  recording  for  deconvolution  pro¬ 
cessing  to  enhance  resolution,  and  wave  migration 
to  transform  distance-time  records  into  distance- 
depth  sections  to  give  a  more  accurate  spatial  image. 

RADAR  EQUIPMENT 
General  operation 

The  radar  used  here  comprises  a  Geophysical 
Survey  Systems,  Inc.,  SIR  Model  4800  control  unit, 
several  different  antennas,  cables,  a  GSSI  DT6000 
digital  tape  recorder  and  a  power  supply.  The  con¬ 
trol  unit  keys  the  transmitter  on  and  off  at  50  kHz 


(synchronized  with  the  receiver),  and  sets  the  scan  known  dielectric  contrasts.  The  scan  shows  a  trigger 

rate  (rate  at  which  echo  scans  are  compiled,  gener-  pulse  (labeled  SOS  for  Start  Of  Scan),  the  Direct 

ally  25.6  scans/s),  scan  time  range  and  the  Time  Coupling  (DC)  between  transmitter  and  receiver  an- 

Range  Gain  (TRG)  to  be  applied  to  the  scans.  The  tennas,  and  reflections  from  subsurface  interfaces, 

transmitter  antenna  radiates  a  broadband  wavelet  The  received  wavelet  forms  in  Figure  1  are  close  to 

that  lasts  from  only  a  few  to  tens  of  nanoseconds  the  transmitted  form  and  last  only  about  two  and 

(ns).  A  separate,  identical  receiver  antenna  is  em-  one-half  cycles.  The  instantaneous  frequency,  de¬ 
ployed  because  echoes  can  return  from  near-surface  fined  as  the  inverse  of  the  time  duration  of  one  cycle, 

events  before  the  transmitter  antenna  has  stopped  is  about  375  MHz  for  the  Model  3102  and  close  to 

radiating.  The  received  signals  are  converted  by  100  MHz  for  the  Model  3107.  For  these  studies 

sampling  into  an  audio  frequency  facsimile  for  fil-  where  the  signal  levels  were  well  above  noise,  the 

tering,  amplifying  and  recording.  phase  of  an  event  is  defined  as  the  sequence  of 

phase  polarities  for  the  first  three  half-cycles  (la- 
Antennas  beled  +  or  -  in  Fig.  1);  other  studies  have  used  only 

All  antennas  are  resistively  loaded  dipoles  that  the  more  visible,  major  half  cycles  (Arcone  et  al.,  in 

radiate  a  wavelet  lasting  about  2  to  2.5  cycles.  The  prep.1).  The  contrast  in  dielectric  constant  e  across 

400-MHz  antenna  pair  are  flared  dipoles  with  a  ta-  the  reflecting  interfaces  determines  whether  the  se- 

pered  resistive  loading,  backshielded  and  housed  in  quence  is  H - 1-  or  — I — .  Reflections  from  interfaces 

one  unit  (Model  3102,  GSSI,  Inc.)  at  a  separation  of  of  known  dielectric  contrasts,  as  well  as  the  DC 

about  15  cm.  The  400-MHz  transmitter  is  rated  by  wavelet  form,  can  be  used  as  a  phase  reference 

the  manufacturer  at  8-W  peak  power  excitation.  The  (Arcone  et  al.,  in  prep.1).  A  wavelet  reflected  from 

100-MHz  antenna  pair  (Model  3107,  GSSI,  Inc.)  are  an  interface  containing  a  low-over-high  e  structure 

of  the  same  design,  but  each  antenna  is  housed  in  a  (e.g.,  ice  over  water)  has  the  same  H - (-  phase  se- 

separate  backshielded  unit.  The  units  were  electri-  quence  as  the  DC  and  that  from  a  high-over-low 

cally  shorted  together,  with  the  antenna  centers  (e.g.,  ice  over  air)  has  a  -  +  -  sequence.  Relative 

separated  by  1  m.  wavelet  phase  will  be  used  to  determine  the  nature 

of  some  of  the  reflectors  within  the  ice  and  snow. 

Waveforms  and  phase  polarity  In  view  of  the  wavelet  duration  and  the  need  for 

A  scan  (Fig.  1)  from  a  400-MHz  profile  presented  separation  of  top  and  bottom  reflections,  resolution 

later  shows  events  received  from  interfaces  between  of  a  thin  air  or  water  layer  embedded  in  ice  would 

require  thicknesses  of  75  and  8  cm  respec¬ 
tively  (eajr  =  1,  Ewater^  87  at  0°C)  for  the  400- 
MHz  antennas,  and  300  and  32  cm  at  100 
MHz.  Practically,  far  thinner  layers  can  be 
resolved  in  a  radar  profile  by  observing 
phase  continuity  from  scan  to  scan  as  the 
overlapping  bottom  reflection  changes 
depth  with  distance.  Resonance  from  layer¬ 
ing  can  be  characterized  by  the  attachment 
of  extra  half-cycles  to  the  primary  layer  re¬ 
flection.  The  attachment  of  an  extra  half¬ 
cycle  at  400  MHz  would  correspond  to  a  19- 
cm  air  layer  and  a  2-cm  water  layer.  Many  of 
the  reflections  seen  later  have  no  extra  oscil¬ 
lations  (i.e.,  resonance)  as  would  occur  in  the 
response  to  an  air  layer  of  less  than  10-cm 
thickness. 

Figure  1.  GPR  400-MHz  wavelets  reflected  from  a  low  density  in -  Antenna  directivity 

elusion  and  brine  within  ice.  Event  SOS  is  an  artificial  start  of  scan  Theoretical  transmitted  radiation  directiv- 
signal;  the  ice  surface  is  at  0  ns.  100-MHz  wavelets  have  similar  form.  ity  f°r  horizontal  point  dipoles  on  the  sur- 

The  sequences  of  half-cycle  phase  polarities  (-\ - 1-  or  -  +  -)  indicate  the  face  a  dielectric  half-space  have  been  dis- 

dielectric  contrast  across  the  reflecting  interface.  The  scan  has  been  high-  cussed  by  Annan  (1973),  Annan  et  al.  (1975), 
pass  filtered.  Engheta  et  al.  (1982)  and  Smith  (1984).  This 
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Figure  2.  Monostatic  radar  transmit-receive  two-way 
directivity  pattern  for  a  finite  size  dipole  operating 
over  ice  or  snow  of  different  dielectric  permittivities 
and  equivalent  densities  p.  The  3-dB  angular  width  of 
the  directivity  in  the  E-plane,  which  is  perpendicular  to 
the  direction  of  surveying,  gives  an  approximate  base  for 
calculating  the  antenna  width  of  sensitivity. 


theory  shows  the  directivity  to  be  nonuniform,  with 
nulls  and  lobes  determined  by  the  value  of  the  criti¬ 
cal  angle  \(/c  =  ±  sin-1  (1  /  n),  measured  from  vertical, 
where  n=Ve  and  is  the  medium  refractive  index.  Di¬ 
rectivity  nulls  in  the  plane  parallel  to  the  antenna 
axis  (electric  field  plane,  or  E-plane),  and  peak  val¬ 
ues  in  the  plane  perpendicular  to  the  antenna  axis 


(magnetic  field  plane,  or  H-plane)  occur  at  \|/c 
which  for  solid  ice  ( n  =  1.78)  is  ±  34°. 

The  transmit-receive  directivities  for  a  horizon¬ 
tal,  unshielded,  finite-size  radar  antenna  operating 
monostatically  (i.e.,  two  antennas  with  zero  offset) 
on  the  surface  of  a  homogeneous  half-space  with 
different  values  of  e  are  shown  in  Figure  2.  The  E- 
plane  is  perpendicular  to  our  profile  directions  and 
the  H-plane  is  parallel.  The  directivities  were  com¬ 
puted  for  the  far-field  electric  field  strength  by  using 
the  superposition  of  fields  radiated  by  discrete  an¬ 
tenna  elements  (Arcone  1994)  and  result  from  com¬ 
puting  either  the  total  energy  content  of  the  pulse 
(area  under  the  waveform)  or  its  peak  value;  both 
give  nearly  identical  results.  The  patterns  model  the 
radiation  response  of  the  Model  3102, 400-MHz  an¬ 
tennas  to  an  idealized,  2.5-ns  (400-MHz)  half-cycle 
current  pulse  that  progressively  attenuates  along 
each  11-cm  half-length  of  the  antenna;  directivity 
patterns  for  the  scaled,  larger  size  100-MHz  anten¬ 
nas  are  nearly  identical. 

The  patterns  of  Figure  2  are  more  directive  than 
for  a  single  transmitter  antenna  owing  to  multiplica¬ 
tion  of  the  identical  transmit  and  receive  directiv¬ 
ities,  and  may  be  considered  the  radar  response  to 
an  isotropic,  polarization-insensitive,  nondispersive 
point  target.  The  two-way  directivity  in  the  E-plane 
shows  a  3-dB  beamwidth  that  progressively  widens 
and  sidelobes  that  progressively  strengthen  as  e  in¬ 
creases.  The  unorthodox  H-plane  pattern  has  maxi¬ 
mum  intensity  in  two  scalloped  lobes,  and  an  effec¬ 
tive  beamwidth  is  not  simply  characterized.  The 
single  antenna  gain  in  ice  (e  =  3.2)  is  about  7  dB. 

FIELD  PROCEDURES  AND 
DATA  PROCESSING 

Field  procedures 

We  did  the  road  and  runway  surveys  by  towing 
the  antennas  behind  a  van  (Fig.  3).  At  the  Williams 
Field  support  facilities,  profiling  was  done  by 
slowly  dragging  the  antennas  over  the  surface  by 
hand  along  measured  transects  within  a  grid.  All 
surveys  were  conducted  with  the  antennas  polar¬ 
ized  perpendicular  to  the  transect  direction  so  that 
the  more  intense,  but  unorthodox,  H-plane  pattern 
was  in  the  line  of  the  transect.  Event  markers  were 
placed  artificially  on  the  records  at  measured  inter¬ 
vals,  generally  every  5  m  for  the  shorter  surveys  at 
Williams  Field.  Scan  durations  were  limited  by  the 
depth  of  the  noise  floor  visible  on  the  oscilloscope 
incorporated  in  the  control  unit,  but  were  usually  set 
to  encompass  the  features  of  interest.  We  used  400  to 


3 


Figure  3.  100-MHz  antennas  in  tow  behind  a  van.  No  reflections  from  the  van  were  apparent  in  the  records.  Van  speed 
was  generally  about  2  m/s. 


750  ns  for  the  100-MHz  studies  and  50  to  200  ns  at 
400  MHz.  A  round-trip  propagation  time  of  100  ns 
corresponds  to  8.4-m  depth  in  solid  (bubble-free) 
ice. 

Data  recording 

All  data  were  recorded  at  a  scan  rate  of  25.6/ 
second  and  scan  sampling  densities  of  512  or  1024  8- 
bit  words.  Both  an  overall  gain  and  a  TRG  were  ap¬ 
plied  to  all  scans.  The  TRG  is  used  to  suppress  the 
DC,  amplify  later  returns  and  keep  all  signals  within 
a  level  required  for  efficient  digitization.  The  GSSI 
Model  4800  TRG  function  is  fairly  fixed  in  form  and 
roughly  compensates  for  the  signal  energy  loss  from 
spherical  wavefront  spreading.  The  TRG  function 
does  not  significantly  alter  any  waveforms  in  this 
study  because  their  duration  is  always  a  small  frac¬ 
tion  of  any  scan. 

Data  processing 

All  data  were  compressed,  or  stacked,  from  five- 
to  twenty-four-fold.  Stacking  suppresses  random 
noise  and  marginally  affected  the  spatial  detail  of 
the  more  intensive,  shorter  profiles,  for  which  the 
stacked  horizontal  data  density  was  about  one  scan 
per  10  cm.  The  lengths  of  the  records  between  dis¬ 
tance  markers  were  equalized  (i.e.,  "justified")  to 
compensate  for  changes  in  towing  speed.  High-pass 
frequency  filtering  was  used  to  alleviate  scan-to- 


scan  gain  variations  resulting  from  inconsistent 
antenna-ground  coupling.  Low-pass  filtering  was 
used  to  alleviate  high-frequency  noise.  Many  pro¬ 
files  were  further  processed  using  predictive  de- 
convolution  to  decrease  the  GPR  wavelet  length  and 
thus  increase  resolution.  This  process  does  not  gen¬ 
erally  retain  the  original  phase  structure  of  the 
wavelet.  All  filtering  and  migration  (discussed  be¬ 
low)  used  algorithms  implemented  by  Boucher  and 
Galinovsky  (1990). 

Echo  time  delay  is  generally  transformed  into 
depth  on  the  basis  of  the  simple  echo  delay  formula 
for  flat  interfaces  or  point  reflectors 

d  =  ct/2n  (1) 

where  d  =  depth  of  a  reflector  (cm) 
t  =  echo  time  delay  (ns) 
c  =  speed  of  electromagnetic  waves  in  a 
vacuum  (30  cm/ns). 

The  factor  of  two  accounts  for  the  round-trip  propa¬ 
gation  path.  Many  reflection  profiles  were  migrated 
to  transform  the  distance-time  record  of  the  raw 
data  into  the  more  meaningful  distance-depth 
record  because  of  irregular  interface  geometry  and 
the  presence  of  obscuring  diffractions.  The  migra¬ 
tions  used  a  single-layer  Kirchhoff  summation  (e.g., 
Yilmaz  1987),  an  approach  commonly  used  in  glaci- 
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Figure  4.  Wiggle  trace  display  of  a  received  scan  and  its  equivalent 
line  intensity  display  should  the  scan  remain  unchanged  over  a 
short  distance. 
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Figure  5.  Approximate  location  of  Williams  Field  and  the  Pegasus 
runway  (adapted  from  Blaisdell  et  al.  1992  and  Mellor  1993).  Several 
sites  were  investigated  at  the  airfield.  Indicated  boundaries  for  the  ice  shelf 
and  transition  zone  are  approximate  also. 
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ology  for  single  interface  migration.  Use  of 
the  single-layer  approach  is  justified  for  sub¬ 
surface  conditions  that  are  approximately  ho¬ 
mogeneous  matrices  of  ice  with  local  inclu¬ 
sions.  The  approach  will  not  work  if  the  in¬ 
clusions  are  so  dense  that  n  is  significantly 
elevated  for  a  large  portion  of  the  ice.  Fea¬ 
tures  of  the  migration  algorithm,  the  effect  of 
the  antenna  directivity  upon  migration,  and 
the  use  of  migration  for  resolving  subsurface 
glacial  features  are  discussed  by  Arcone  et  al. 
(in  prep.1) 

Successful  migration  of  a  single  interface 
depends  on  accurate  knowledge  of  the  veloc¬ 
ity,  or  n,  of  the  overlying  layer.  In  the  absence 
of  drilling  to  obtain  true  reflector  depth, 
single  layer  velocity  can  be  determined  at  the 
surface  either  from  surface  refractions  prop¬ 
agated  between  continuously  separated 
antennas  (Arcone  1984),  from  the  asymp¬ 
totic  slopes  of  diffractions  originating 
within  or  at  the  bottom  of  the  layer,  or  by 
using  a  dual  antenna  arrangement  to  solve 
for  both  depth  and  n  simultaneously 
(Kovacs  et  al.  1982).  The  first  two  methods 
were  used  here,  but  only  the  latter  two 
give  reasonably  accurate  results  when  n 
continually  increases  with  depth,  as  is  the 
case  for  snow,  because  surface  refract¬ 
ions  propagate  mainly  in  the  near-surface 
layers.  Deep  reflections  and  diffractions 
will  spend  a  disproportionate  amount  of 
time  within  the  upper  layers  because  of 
the  gradient  of  n,  and  so  will  give  an  ef¬ 
fective  n  biased  toward  the  surface.  Mi¬ 
grations  in  these  studies  used  a  value  nm 
calculated  from  matching  the  slopes  of  the 
diffractions'  asymptotes,  generally  to  an 
accuracy  of  about  ±2%.  This  procedure 
produced  good  migrations  of  the  diffrac¬ 
tions,  but  nm  is  not  necessarily  the  aver¬ 
age  refractive  index  of  the  section  mi¬ 
grated. 

Data  are  displayed  as  echo-time-of- 
retum  vs.  profile  distance  using  either 
wiggle  trace  or  line  intensity  format  as  de¬ 
picted  in  Figure  4.  A  nonlinear  scheme  of 
scale  intensity  that  ■  emphasizes  weaker 
events  is  used  for  the  line  intensity  format. 
Profiles  were  printed  on  an  Hewlett-Pack¬ 
ard  Paintjet  printer.  The  quality  of  the 
printout  was  generally  very  good,  but 
inferior  to  the  display  on  the  computer 
monitor. 
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SITE  LOCATION 

All  sites  investigated  were  on  the  Ross  Ice  Shelf 
near  McMurdo  Station,  Antarctica,  and  were  associ¬ 
ated  with  operations  at  Williams  Field  and  the  Pega¬ 
sus  runway  (Fig.  5).  Profiles  were  recorded  along  the 
Pegasus  runway,  the  access  road  to  the  Pegasus  run¬ 
way,  and  three  sites  at  the  Williams  Field  support 
facilities.  All  sites  were  located  in  topographically 
depressed  areas  that  were  formed  by  maintenance 
activities  that  have  removed  snow,  some  continu¬ 
ously  for  several  years.  The  surveys  were  performed 
between  14  and  18  January  1993.  Weather  during 
profiling  was  partially  cloudy  to  overcast,  with  light 
winds  and  temperatures  around  0°C.  Surface  melt¬ 
ing  is  common  on  sunny  days  at  the  peak  of  the 
austral  summer  and  small  meltwater  pools  were 
encountered  on  the  Pegasus  runway  and  occasional 
wet  snow  at  the  Williams  Field  solid  waste  dump 
site. 

Both  the  Pegasus  runway  (probably)  and  the 
Williams  Field  area  are  within  the  accumulation 
zone  where  snow  is  accumulating  on  top  of  glacial 
ice  (Fig.  5).  The  Pegasus  site,  however,  is  in  a  transi¬ 
tion  zone  of  superimposed  ice  that  is  refrozen  water 
from  runoff  and  melted  snow  (Blaisdell  et  al.  1992), 
whereas  Williams  Field  is  in  a  zone  of  permanent 
snow.  Cores  taken  by  Blaisdell  et  al.  along  the  Pe¬ 
gasus  runway  to  2.5-m  depth  reveal  bubbly  layers 
with  densities  as  low  as  0.6  g/cm3.  Subsurface  melt 
pools  1  to  1.5  m  deep  have  been  observed  farther  to 
the  west  within  glacial  "blue"  ice  and  are  attribut¬ 
able  to  the  greenhouse  effect  (Paige  1968).  The  same 
process  is  known  to  form  near-surface  melt  pools  at 
Pegasus.  They  are  generally  prevented  by  roughen¬ 
ing  the  surface.  The  accumulation  zone  has  been 
well  investigated,  but  particularly  by  Kovacs  and 
Gow  (1975)  and  Kovacs  et  al.  (1982),  who  used  radar 
and  coring  to  study  the  infiltration  of  seawater  into 
the  snow;  further  discussion  of  their  results  is  given 
later.  We  assume  that  snow  exists  at  depth  beneath 
the  Pegasus  site,  based  on  the  continuity  of  radar 
profiles  of  brine  infiltration  given  below,  but  we  do 
not  have  deep  enough  cores  from  this  area  to  con¬ 
firm  this.  Drilling  in  December  1991  by  the  Polar  Ice 
Coring  Office  (PICO)  revealed  the  ice  shelf  to  be  ap¬ 
proximately  33  m  thick  at  the  north  end  of  Pegasus 
runway.  Freeboard  above  sea  level  was  4.8  m. 


RESULTS 
Pegasus  runway 

The  objective  of  this  survey  was  to  determine  if 
GPR  could  detect  voids  or  water  pockets  near  the 
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surface,  as  these  could  cause  failure  of  the  ice  pave¬ 
ment  during  wheeled  landings  of  heavy  cargo  air¬ 
craft.  Figure  6  depicts  the  runway  and  the  place¬ 
ment  of  the  profiles,  all  of  which  were  run  on  the 
west  side.  The  runway  was  flagged  every  500  ft 
(152.4  m)  over  a  distance  of  10,500  ft  (3200  m).  The 
100-MHz  transducer  unit  was  first  used  with  a  time 
range  set  to  a  depth  of  58  m  (for  solid  ice;  greater  for 
snow).  The  profile  revealed  only  a  few  diffractions 
originating  later  than  a  strong  horizon  at  about  9-m 
depth.  Subsequently,  400-MPlz  profiles  were  run 
and  only  they  are  presented  as  they  give  much 
greater  detail. 

Figure  7  gives  a  compressed  profile  of  the  Pe¬ 
gasus  runway  from  11,000  to  0  ft  (3353  to  0  m).  At  all 
depths  and  distances,  the  slopes  of  the  diffraction 
asymptotes  give  an  e  very  near  3.2,  the  value  for 
solid  (bubble-free)  ice,  and  which  agrees  with  the  ice 
core  observations  of  Blaisdell  et  al.  (1992)  of  pre¬ 
dominantly  solid  ice  in  the  upper  2  m.  The  domi¬ 
nant  feature  is  the  reflections  originating  between 
about  8-  and  9-m  depth.  We  believe  that  this  radar 
horizon  is  brine  infiltrating  permeable  snow  at 
depth  because  it  rises  toward  the  north  end  (where 
the  shelf  thins)  and  continues  into  the  accumulation 
zone  (next  section),  where  drilling  has  confirmed 
this  interpretation  (Kovacs  et  al.  1982).  In  tire  Pe¬ 
gasus  area,  the  near-surface  material  is  classified  as 
superimposed  ice,  which  is  refrozen  meltwater  from 
snow  and  runoff.  Alternatively,  the  horizon  could  be 
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Figure  7.  Compressed  profile  at  400  MHz  of  the  west  side  of  the  Pegasus  runway. 


brine  infiltration  within  permeable  layers  of  the  su¬ 
perimposed  ice,  or  accreted  sea  ice/  which,  there¬ 
fore,  must  be  more  than  20  m  thick.  A  deep  core 
would  resolve  this  issue,  but  this  was  not  the  pur¬ 
pose  of  our  survey. 

Above  the  8-  to  9-m  horizon  in  Figure  7  is  a  more 
diffuse  level  of  reflections,  with  many  horizons  con¬ 
centrated  mainly  between  3  and  5  m  depth.  These 
reflections  are  caused  by  density  inhomogeneities  in 
the  ice,  as  discussed  below.  The  time  delay  between 
the  3-  to  5-  and  8-  to  9-m  levels  precludes  interpreta¬ 
tion  of  the  8-  to  9-m  reflections  as  multiples  from  the 
upper  horizons.  Reflections  from  the  two  levels  are 
also  distinctly  different  in  structure,  and  in  phase 
(discussed  next). 

More  detailed,  uncompressed  profiles  of  runway 
sections  are  shown  in  Figure  8.  These  section  pro¬ 
files  (as  does  the  rest  of  the  whole  profile)  reveal 
distinct  linear  horizons  up  to  30  m  long  between  the 
3-  and  7-m  depths.  They  are  probably  responses  to 
bubble  layers  incorporated  during  melting  and  re¬ 
freezing.  Even  greater  resolution  of  the  upper  level 
events  is  seen  in  the  100-ns  time  range  profile  shown 
in  Figure  9.  At  this  range  many  events  appear  as 
single  surface  reflections,  as  would  occur  for  dielec¬ 
trically  thin  density  anomalies.  The  brine  "reflec- 


*  Personal  communication  with  Dr.  A.J.  Gow,  CRREL, 
1994. 


tions"  at  8-9  m  in  Figure  8  are  actually  composed  of 
diffractions  from  discontinuities  in  the  dielectric 
structure.  These  diffractions  imply  that  the  brine  is 
intruding  along  discrete  pathways. 

As  discussed  earlier,  reflection  phase  reversals 
relative  to  the  direct  coupling  between  these  closely 
spaced  transmit  and  receive  antennas  are  ascribable 
to  decreases  in  e  with  depth  across  a  reflecting  inter¬ 
face.  Many  of  the  upper  level  horizons  in  Figure  8 

generally  have  an - t-  structure  (white-black- white 

banding  in  the  profile — more  apparent  when  the 
records  are  viewed  edge  on)  to  the  strongest  and 
most  visible  half-cycles  of  the  wavelet,  which  indi¬ 
cates  a  negative  dielectric  anomaly  (lower  density). 
In  contrast,  the  horizon  at  about  the  9-m  depth  has  a 
— i —  wavelet  structure  (Fig.  1),  which  indicates  a 
positive  dielectric  anomaly.  This  interpretation  is 
verified  by  a  5-m  ice  core  (Fig.  10)  obtained  at  the 
10,492-ft  (3198-m)  location  on  the  runway  where  re¬ 
flections  occur  in  Figure  9  at  about  the  1.9-  and  3.5-m 
depths.  The  core  density  anomalies  are  associated 
with  air  bubbles  and  occur  at  the  depths  of  reflec¬ 
tions  in  an  individual  scan  from  this  station,  also 
shown  in  Figure  10.  Events  1  and  2  in  the  scan  have 
phase  sequences  opposite  that  of  the  direct  cou¬ 
pling.  Event  3  has  the  same  phase  sequence  and 
may  be  a  secondary  reflection  from  the  bottom  of 
the  layer  at  3.5  m.  A  secondary  reflection  from  this 
lower  layer  would  be  visible  because  the  layer  is 
thicker  than  the  layer  at  1.9  m,  as  seen  in  the  core 
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Figure  8.  Uncompressed  profiles  at  400  MHz  of  two  500ft  (152-m)  segments  along  the  west  side  of 
Pegasus  runway's  south  end.  The  horizons  beginning  at  about  3  m  depth  appear  to  be  bubble  layers.  The 
phase  sequence  of  the  major  half- cycles  of  the  wavelets  is  generally  +  -  +  (white-black-white  banding)  and 
is  opposite  to  that  of  the  brine  reflections  (—  +  —)  at  9  m  depth. 
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Figure  9.  Uncompressed  profile  at 
400  MHz  showing  the  density  an¬ 
omalies  near  the  10,500-ft  (3200-m) 
distance.  Arrows  point  to  the  anoma¬ 
lies  investigated  in  Figure  10.  Some  of 
the  internal  reflections  persist  to  30  m 
length. 


8 


P  e 

3 

Density  (g/cm  )  Relative  Dielectric  Permittivity 


Signal  Amplitude 


Figure  10.  Computed  density  and  e  profiles  of  a  5-m  core,  and  a  radar  profile  scan  taken  at 
the  10,492-ft  (3198-m)  distance  along  Pegasus  runway.  Density  is  determined  by  weighing 
and  measuring  sample  volume  and  is  accurate  to  about  2%.  Values  of  e  based  on  the  volumetric 
ratios  of  ice  and  air  (e.g.,  Gumming  1952).  Density  is  generally  above  0.85,  verifying  the  near¬ 
solid  ice  interpretation  ofe  from  the  diffractions.  The  phase  of  the  reflected  wavelets  relative  to  the 
direct  coupling  indicates  a  transition  from  high  to  low  e  at  the  reflecting  interfaces  (see  Fig.  1 
also). 
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Figure  11.  Uncom¬ 
pressed  profile  along 
Pegasus  runway.  Near¬ 
surface  melt  pools  occur 
between  8000  and  8100 
ft  (2.438  and  2.469  km). 
The  record  has  been  hori¬ 
zontally  filtered  to  elimi¬ 
nate  the  surface  coupling 
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profile.  The  separation  of  events  2  and  3  is  4  ns, 
which  corresponds  to  a  layer  thickness  of  about  40 
cm,  given  an  e  of  2.1,  and  is  consistent  with  the 
anomaly  width  in  the  density  profile. 

Of  particular  concern  to  Pegasus  operations  are 
voids  or  melt  pools  within  about  1  m  of  the  surface. 
Several  strong  and  extensive  near-surface  reflections 
were  seen  along  the  Pegasus  profile  between  8000 
and  10,000  ft  (2438  and  3048  m),  some  of  which  are 
shown  in  Figure  11.  These  are  meltwater  pools  and 
they  give  a  characteristic,  diffraction-dominated  re¬ 
sponse  for  localized,  near-surface  moisture.  From  0 


to  8000  ft  (0  to  2438  m),  there  was  near-surface  melt¬ 
water  along  the  profile  at  158  locations,  but  each  was 
limited  to  puddles  less  than  about  20  cm  in  extent. 
Therefore,  about  1%  of  the  profile  encountered  near¬ 
surface  or  surface  meltwater.  An  airborne  or  plat- 
form-mounted  microwave  survey  would  be  better 
suited  to  this  application,  as  all  indications  of  melt¬ 
water  appear  to  originate  very  near  or  at  the  surface. 
The  existence  of  these  small  but  weak  areas  in  the  ice 
of  the  runway  is  reason  to  continue  to  proof-test  the 
ice  pavement  periodically  with  a  load  cart. 
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Figure  12.  Profile  sections  (PEG-A-B-C-D-WIL)  along  the  access  road  from 
Williams  Field  to  the  north  end  of  Pegasus  runway.  Line  marked  Flight  94  is  the 
approximate  profile  done  by  air  in  January  1994  and  discussed  in  Appendix  A. 


Table  1.  Ice  shelf  thicknesses  at 
locations  shown  in  Figure  12. 

Location  Thickness  (m) 

Pegasus  North 
A 
B 
C 
D 

H  (Williams  Field) 

Gow  (1975)  in  using  the  ex¬ 
pression  brine  "layer/7  as  they 
point  out  that  it  cannot  pen¬ 
etrate  snow  with  a  density 
greater  than  about  0.8.  How¬ 
ever,  we  can  also  call  it  a  zone 
of  brine  permeation  or  pen¬ 
etration. 

The  100-MHz  antenna  was 
used  because  the  brine  layer 
eventually  deepens  to  the  point 
where  sample  spacing  within 
the  radar  scans  would  have 
been  insufficient  to  reproduce 
the  reflection  waveforms  accu¬ 
rately  at  400  MHz.  The  accuracy 
of  the  distances  shown  for  this 
and  all  other  sections  is  esti¬ 
mated  at  about  ±20  m. 
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Access  road  from  Williams  Field 
to  Pegasus  runway 

The  purpose  of  this  survey  was  to  profile  the 
depth  to  the  brine  layer  and  to  examine  the  access 
road  for  any  unusual  features,  such  as  cracking  in 
the  ice.  The  approximate  routing  of  the  access  road 
is  shown  in  Figure  12.  The  road  is  divided  into  sec¬ 
tions,  the  ends  of  which  are  indicated  by  letters.  Re¬ 
cent  hot  water  drilling  (December  1991)  by  PICO 
has  determined  the  ice  shelf  thicknesses  presented 
in  Table  1. 

A  compressed,  100-MHz  profile  for  the  access 
road  section,  progressing  from  the  north  end  of  the 
Pegasus  runway  to  point  A  in  Figure  12,  is  shown  in 
Figure  13.  The  strong  reflection  seen  at  8-9  m  depth 
beneath  the  Pegasus  runway  continuously  extends 
through  this  record  and  on  to  Williams  Field,  an  area 
in  which  the  horizon  is  known  to  be  intruding  sea¬ 
water  as  discussed  above.  Only  a  vertical  portion  of 
the  access  road  records  is  shown  because  there  are 
no  other  primary  reflections  from  beneath  the  brine 
layer,  brine  being  highly  absorbing  of  radio  wave 
energy.  We  follow  the  terminology  of  Kovacs  and 


Although  it  is  not  apparent  in 
the  compressed  profiles,  the  brine  layer  reflections 
themselves  are  not  actually  continuous  reflections, 
but  composites  of  multiple  diffractions,  as  was 
seen  at  the  Pegasus  runway.  These  diffractions  are 
scatter  from  discrete  areas  and  show  that  the  brine 
layer  seems  to  intrude  along  discrete  pathways 
within  the  snow.  No  diffractions  originating  at  or 
above  the  brine  layer  that  might  be  characteristic  of 
cracking  or  crevassing  are  apparent  in  any  of  the 
access  road  profiles. 

The  depth  scale  at  the  start  of  Figure  13  is  based 
on  nm=1.78  (e  =  3.2),  which  is  determined  by  slopes 
of  the  diffraction  asymptotes  at  the  north  end  of  the 
Pegasus  runway  where  the  profile  starts.  This  solid 
ice  value  of  e  gives  the  minimum  possible  depth  to 
the  brine  layer  and  is  consistent  with  the  density  of 
ice  cores  taken  near  Pegasus  (Blaisdell  et  al.  1992). 
The  brine  level  depth  varies  from  about  8  to  nearly 
10  m  and  stays  at  this  level  just  beyond  point  C 
along  the  road,  where  it  will  be  seen  later  that  this 
level  is  being  superseded  by  a  younger  and  more 
shallow  brine  intrusion.  Figure  13  also  shows  a  dif¬ 
ferent  horizon  dipping  from  southwest  to  northeast. 
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Figure  13.  Twenty-five-fold  compressed,  100-MHz  radar  profile  from  the  north  end  of  Pegasus  runway 
to  point  A.  Horizontal  noise  bands  have  been  removed.  The  dominant  feature  is  the  returns  from  the  brine 
layer.  The  origin  of  the  diffuse  horizon  descending  to  the  northeast  in  the  first  half  of  the  profile  is  not  known, 
but  it  cannot  be  the  contact  between  snow  and  ice  because  there  is  no  surface  snow  at  Pegasus. 
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Figure  14.  Deconvolved,  100-MHz  profile  of  sections  A-B  and  B-C. 


The  origin  of  this  is  not  known,  but  it  cannot  be  a 
snow-over-ice  contact  because  there  is  no  near-sur- 
face  snow  at  Pegasus. 

A  deconvolved,  100-MHz  profile  for  sections  A- 
B  and  B-C  of  the  road  is  shown  in  Figure  14.  The 
profile  contains  two  features  of  interest:  folded  stra¬ 
tification,  and  two  apparent  brine  layers  with  dis¬ 
continuities  between  A  and  B.  The  depth  scale  is 
based  on  diffraction  slopes  at  point  B  that  corre¬ 
spond  to  nm=1.78  (e  =  3.2),  but  there  is  evidence  to 
show  dielectric  variability.  Refraction  profiles  be¬ 
tween  A  and  B  give  a  surface  value  n=1.4,  which 
corresponds  to  a  density  of  0.5.  At  point  B,  the  pen¬ 
etration  rate  of  a  hot  point  probe  decreased  (because 
of  heat  absorption  by  the  brine)  rapidly  at  9.3  m; 
probing  between  the  discontinuities  is  discussed  be¬ 


low.  The  radar  profile  gives  a  depth  to  brine  of  10.0 
m  at  point  B.  Using  the  hot  point  data  for  calibration 
suggests  that  n=1.9,  which  is  slightly  high.  The  dis¬ 
crepancy  cannot  be  explained  without  knowing  the 
exact  gradient  of  the  brine  concentration,  and  the 
exact  ray  path  through  the  dielectric  stratification. 
The  folding  suggests  that  ray  paths  to  the  brine  layer 
were  not  always  direct,  but  the  correction,  discussed 
next,  may  be  minor. 

A  deconvolved  and  expanded  100-MHz  profile 
of  the  folding  of  density  stratification  in  section  B-C 
is  shown  in  Figure  15.  Scale  compression  and  a  hori¬ 
zontal  (or  "background  removal")  filter  to  remove 
the  resonance  of  the  direct  coupling  have  been  used 
to  make  the  folding  more  visible.  The  stratification 
shows  that  the  gradient  of  n  is  discontinuous;  most 


11 


Figure  15.  Expanded  section  of 
the  profile  in  Figure  14,  show¬ 
ing  detail  of  folding  in  section 
B-C.  Background  noise  has  been 
removed  to  show  the  diffractions 
that  constitute  the  brine  level  dis¬ 
continuity  at  10  m  depth. 
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Figure  16.  Compressed,  deconvolved  and  migrated  profile  showing  detail  of  the  brine 
discontinuities  and  an  exposure  of  a  lower  brine  layer  seen  in  section  A-B  in  Figure  14. 
The  migration  used  e  =  3.2.  The  diffuse  layer  at  about  the  10-m  depth  probably  contains  brine. 


1.8 


of  the  change  must  occur  near  the  surface  where 
n  -  1.4  because  nm  is  between  1.70  and  1.78.  The 
width  of  one  limb  of  a  fold  is  about  150  m.  This  is 
five  times  the  brine  depth  and  corresponds  to  an 
angle  of  only  about  4°  from  horizontal.  This  slope, 
along  with  the  gradual  change  in  n  with  depth, 
forces  ray  path  bending  to  provide  normal  incidence 
at  the  brine  layer  and  reflection  back  to  the  receiver 
antenna.  Hie  total  bending  is  less  than  2°,  provides 
less  than  1%  correction  to  the  ray  path  length  and 
does  not  affect  nm  significantly. 

A  deconvolved  migration  of  the  large  disconti¬ 
nuity  containing  two  diffuse  radar  horizons  within 
section  A-B  is  shown  in  more  detail  in  Figure  16. 
Within  the  discontinuity,  returns  from  the  deeper 
horizon  are  stronger  than  those  from  higher  in  the 
snow  so  that  these  deeper  events  cannot  be  multiple 
reflections.  The  migration  used  an  nm  of  1.78  based 
on  the  diffraction  asymptotes,  but  slightly  lesser  val¬ 
ues  work  nearly  as  well  because  precompression, 
needed  to  obtain  an  effective  migration  aperture, 


mitigates  migration  noise.  The  profile  reveals  very 
little  continuity  in  the  upper  horizon,  but  fairly  good 
continuity  in  the  lower.  The  penetration  rate  of  a  hot 
point  probe  half  way  between  the  ends  of  the  dis¬ 
continuity  decreased  rapidly  at  9.6  m  (the  deeper 
layer  was  not  reached),  while  the  radar  record 
shows  this  diffuse  upper  horizon  to  begin  at  about 
9.8-m  depth.  Only  the  southwest  side  of  the  discon¬ 
tinuity  shows  any  transition  between  upper  and 
lower  horizons,  and  both  sides  of  the  upper  horizon 
have  abrupt  terminations.  The  upper  horizon  may 
represent  a  dense  layer  of  snow  into  which  brine 
from  either  side  is  slowly  diffusing,  while  the  lower 
horizon  is  either  an  exposed  section  of  an  older 
brine  layer,  sea  water  or  accreted  sea  ice  at  the  bot¬ 
tom  of  the  shelf.  Ice  thicknesses  measured  a  year 
earlier  were  40  m  at  point  A  and  36  m  at  point  B, 
which  suggests  that  the  lower  horizon  is  not  the  bot¬ 
tom  of  the  shelf. 

Kovacs  and  Gow  (1975)  and  Kovacs  et  al.  (1982) 
give  a  thorough  discussion  of  ice  shelf  breakout. 
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Figure  17.  100-MHz  deconvolved  profile  of  sections  C-D  and  D-WIL.  The  depth  scale  is 
based  on  e  =  2.7. 
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Figure  18.  Time  (bottom)  and  migrated  profiles  containing  the  brine  discontinuities  seen 
in  sections  C-D  (left)  and  D-WIL  (right).  The  expanded  distance  scale  is  25xless  than  that  of 
Figure  17.  Diffractions  indicate  discontinuities.  Apparent  layering  above  the  brine  level  in  the 
migrated  profile  at  right  is  an  artifact  of  the  migration. 


subsequent  brine  "wave"  intrusion  and  overriding 
of  previous  brine  levels  that  existed  before  the 
breakout.  Their  radar  data  profile  stepped  intru¬ 
sions  as  they  extend  away  from  the  edge  of  the  ice 
shelf.  The  example  of  Figure  16  seems  to  illustrate  a 
cross  section  of  this  process  at  a  location  where  the 
intruding  front  is  not  uniform.  As  mentioned  above 
and  as  will  be  shown  next,  the  10-m  brine  level  seen 
in  Figure  16  is  being  overrun  by  a  younger  and  more 
shallow  brine  level. 

Profiles  for  sections  C-D  and  D-WIL  are  shown 
in  Figure  17.  At  Williams  Field  the  ice  shelf  thickens 
because  of  snow  accumulation  and  the  brine  level 


drops  below  20  m.  Diffraction  slopes  at  this  end  cor¬ 
respond  to  an  nm  of  1.64  (e  =  2.7).  Two  brine  level 
changes  are  seen  here  and  these  are  the  same  type  of 
features  discussed  above,  but  are  more  similar  to  the 
examples  of  Kovacs  and  Gow  (1975).  The  brine  layer 
rises  just  beyond  point  C  because  the  edge  of  the  ice 
shelf  (where  it  contacts  the  sea)  is  closest  between 
points  C  and  D  and  then  falls  after  point  D  as  tire  sea 
ice  contact  becomes  more  distant  (see  Fig.  12).  Both 
discontinuities  are  then  part  of  the  same  brine  front 
between  points  C  and  D.  That  brine  front  is  the  re¬ 
sult  of  calving  in  this  area.  An  airborne  profile  of  this 
section,  approximately  from  points  C  to  WIL,  was 
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run  the  following  year  and  is  shown  in  Ap¬ 
pendix  A. 

Time  and  migrated  sections  of  both  brine 
discontinuities  are  shown  in  Figure  18. 
Only  enough  compression  was  used  to  per¬ 
mit  the  migration  aperture  of  61  scans  (the 
maximum  allowable  by  the  software  for  a 
1024  sampling  density)  to  cover  the  diffrac¬ 
tions  adequately,  but  the  scale  is  still  25x 
less  than  that  of  Figure  17.  In  contrast  to  the 
discontinuities  seen  in  Figure  16,  the  brine 
wave  front  between  points  C  and  D  appears 
gradual,  contains  no  significant  diffractions 
(in  the  time  section)  and  was  probably  tra¬ 
versed  obliquely  to  the  brine  front  normal. 
The  front  of  the  intrusion  between  D  and 
WIL,  however,  contains  at  least  two  discon¬ 
tinuities  as  evidenced  by  the  diffractions  in 
the  time  section  (Fig.  17).  The  migration 
shows  that  between  the  two  breaks  there  is 
an  almost  continuous  and  coherent  sloping 
reflection,  which  indicates  a  smooth  brine 
front  propagating  northeast. 
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Figure  19.  Williams  Field  support  facilities  (not  to  scale). 


Williams  Field 

A  simplified  map  of  Williams  Field  is  shown  in 
Figure  19.  The  layout  shows  the  suspected  location 
of  a  (presumed,  surface)  fuel  spill  near  the  fuel  stor¬ 
age  tanks,  a  suspected  buried  solid  waste  dump  and 
old  and  new  sewage  sumps. 

Fuel  spill  site 

Our  objectives  at  the  spill  site  were  to  locate  and 
delineate  any  fuel  that  may  have  not  been  removed 
from  a  previous  Navy  cleanup  operation.  During 
this  cleanup,  much  snow  was  melted,  which  obvi¬ 
ously  refroze  and  has  subsequently  been  buried. 
The  main  concern  was  that  much  of  the  fuel  had 
percolated  downward  and  was  not  retrieved. 

A  fuel  spill  in  snow  can  provide  a  subtle  dielec¬ 
tric  contrast  with  the  surrounding  snow.  Fuel  has  a 
typical  dielectric  constant  e  =  2.2  and,  when  filling 
the  pore  spaces  in  snow  of  about  0.4  density,  would 
produce  an  e  of  about  2.6.  This  would  not  provide 
much  more  contrast  than  would  solid  ice  (e  =  3.2).  If 
the  uncontaminated  snow  already  has  an  e  =  2.8  or 
higher,  then  it  may  not  be  permeable.  This  was  obvi¬ 
ously  not  the  case  because  the  fuel  infiltrated  the 
snowpack  and  was  not  ponded  on  the  surface. 

The  layout  of  four  profiles  taken  over  the  sus¬ 
pected  fuel  spill  area  is  shown  in  Figure  20.  The 
numbered  boxes  indicate  the  numbered  fuel  storage 
tanks.  The  diametric  layout  of  the  65-  to  80-m-long 
profiles  was  centered  at  the  suspected  origin  of  the 
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Figure  20.  Survey  grid  of  the  fuel  spill  site  at  Williams 
Field ;  the  position  cf  the  fuel  tanks  is  for  orientation  only 
and  not  to  scale.  The  dashed  line  indicates  the  suspected 
spill  area  interpreted  from  the  radar  profiles  and  the  cross 
hatching  indicates  a  dipping  horizon  that  leads  to  a  bright 
reflection  at  7  m  depth  and  70  m  distance  along  B-B'.  The 
dots  indicate  particularly  strong  diffractions,  which  are  at  the 
5-  to  8-m  depths. 
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22.  Deconvolved,  400-MHz  profiles  of  the  fuel  spill  area.  The  depth  scale  corresponds  to  am-  2.3. 


spill.  Profiles  were  made  with  the  400-MHz  trans¬ 
ducer  using  a  time  range  of  200  ns.  However,  within 
this  time  range,  events  that  could  possibly  be  a  re¬ 
sponse  to  a  fuel  spill  occurred  within  the  first  80  ns 
throughout  all  profiles  (Fig.  21).  Each  profile  has 
been  deconvolved  and  has  had  the  antenna  direct 
coupling  removed  to  reveal  any  near-surface  dis¬ 
turbance.  Diffraction  asymptotes  show  values  of 
nm=  1.52  (e  =  2.3),  which  implies  that  any  fuel  spill 
was  not  pervasive  and  could  be  detected.  This  value 
of  e  gives  a  depth  of  about  8  m  to  the  records,  which 
should  be  deeper  than  any  excavation  done  here. 
Migrations  were  not  done  because  the  diffractions 
highlight  the  disturbances. 

All  profiles  in  Figure  21  show  linear  horizons, 
presumed  to  be  density  layering  in  the  snow.  How¬ 
ever,  any  of  these  horizons  could  be  fuel  or  ice,  but  it 
is  not  known  if  fuel  can  concentrate  rather  than  dis¬ 
perse  in  snow.  There  are  particularly  strong  re¬ 


flections  at  about  20, 40  and  60  ns,  which  translate  to 
depths  of  about  2,  4  and  6  m  respectively.  All  the 
profiles  contain  areas  of  strong  diffractions  concen¬ 
trated  in  the  central  area,  and  some  occur  in  mul¬ 
tiples  (e.g.,  at  44  m  in  A'-A  and  at  36  m  in  D-D').  In 
addition.  Line  B-B'  shows  a  weak  horizon  that 
dips  from  40-  to  70-m  distance,  at  which  point  it 
terminates  in  a  particularly  strong  reflection.  The 
area  of  the  diffractions  and  the  dipping  horizon 
along  B-B'  may  be  the  area  of  the  fuel  spill,  the  ex¬ 
tent  of  which  is  delineated  in  Figure  20.  By  the  time 
of  this  writing,  drilling  within  the  interpreted  area  of 
disturbance  in  Figure  20  had  not  been  pursued. 

Solid  waste  dump  site 

The  objective  of  this  limited  survey  was  to  verify 
the  suspected  area  of  a  solid  waste  dump  that  was 
used  for  a  number  of  years  until  1990.  Metal  waste 
provides  strong  radar  targets.  The  layout  and  ap- 
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Figure  21  ( cont'd ).  Deconvolved,  400-MHz  profiles  of  the  fuel  spill  area.  The  depth  scale  corresponds  to  an 
e  =  2.3. 


proximate  location  of  a  single,  100-m-long  profile 
across  the  dump  on  the  north  side  of  Williams  Field 
was  shown  in  Figure  19.  The  profile  was  centered  on 
and  perpendicular  to  a  plowed  swath  approximate¬ 
ly  15  m  wide.  This  swath  was  presumed  to  run 
through  the  suspected  northwest-southeast  orien¬ 
tation  of  the  dump.  The  location  of  the  line  was  not 
referenced  exactly  relative  to  any  fixed  landmarks, 
as  this  was  only  a  test  to  see  if  GPR  could  find  the 
buried  waste,  to  be  followed  up  with  backhoe  exca¬ 
vation  to  prove  the  value  of  the  radar  system  to 
maintenance  personnel  at  McMurdo  Station. 

The  first  100  ns  of  the  400-MHz  profile  of  the 
solid  waste  dump  line  is  shown  in  Figure  22;  no  sig¬ 
nificant  information  was  obtained  at  later  return 
times.  The  upper  part  of  the  record  shows  subtle 
horizons  that  rise  toward  the  center  of  the  record. 
These  horizons  are  assumed  to  be  snow  layers  be¬ 
cause  they  rise  to  where  plowing  had  removed  the 


snow  from  the  central  part  of  the  transect.  Several 
anomalous  returns  appear  beneath  these  horizons, 
with  the  strongest  centered  at  the  20-,  55-  and  87-m 
distances.  All  diffraction  asymptotes  associated 
with  these  disturbances  give  an  nm=1.38  (e  =  1.9),  a 
value  that  corresponds  with  a  snow  density  of  ap¬ 
proximately  0.5.  This  value  of  e  is  most  likely  a  mini¬ 
mum  for  this  material  and  shows  the  influence  of 
the  near-surface  layers  in  refracting  subsurface  dif¬ 
fractions.  The  value  is  used  to  derive  the  depth  scale 
of  the  profile,  which  provides  a  maximum  range  be¬ 
cause  the  snow  density  increases  with  depth. 

Subsequent  to  profiling,  a  section  near  the  55-m 
mark  was  excavated  by  bulldozer  and  metal  debris 
was  found  at  a  depth  of  3.4  m.  This  compares  well 
with  the  bright  reflection  seen  at  about  3.7  m  depth 
in  Figure  22.  Hard  ice  was  encountered  at  the  1.5-m 
depth,  and  explains  the  earlier  radar  horizon  at 
about  1.5  to  2  m  depth  near  the  55-m  mark. 
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Figure  22.  400-MHz  profile  of  the  solid  waste  dump  survey  line.  Upper  horizons  are  presumably 
snow  layers.  The  diffraction  slopes  give  an  e  =  1.9,  a  value  that  gives  a  maximum  depth  range  for  the 
given  time  range.  Excavation  about  the  55-m  mark  revealed  hard  ice  at  about  the  1.5-m  and  metal  debris  at 
about  the  3.4-m  depths. 


Sewage  sumps 

The  objective  of  this  last  survey  was  to  locate  the 
upper  surface  of  sewage-soaked  snow  at  the  buried 
sewage  sumps  (Fig.  19)  and  delineate  its  lateral  ex¬ 
tent.  Snow  soaked  with  liquid  sewage  provides  a 
strong  dielectric  contrast  with  uncontaminated 
snow  because  the  value  of  e  for  water  at  the  freezing 
point  is  about  87,  although  it  is  slightly  lower  when 
filled  with  organic  waste.  The  layout  of  one  longi¬ 
tudinal  and  two  cross-sectional  profiles  taken  at  the 
sewage  sumps  is  shown  in  Figure  23.  The  sumps  are 
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Figure  23.  Plan  view  showing  survey  grid  and  location  of  the 
sewage  sumps  at  Williams  Field.  The  dashed  lines  indicate  the 
zone  of  contaminated  snow  as  interpreted  from  the  radar  data.  Pj 
and  ?2  denote  holes  melted  with  a  hot  point  drill  to  verify  GPR 
findings. 


33  m  apart  in  a  west-southwest-east-northeast  ori¬ 
entation.  The  eastern  sump  is  older,  filled  with  sew¬ 
age  and  no  longer  in  use.  Some  of  the  sewage  there 
has  probably  refrozen.  As  more  and  more  sewage  is 
discharged  into  the  snow,  the  contaminated  zone 
rises.  The  western  sump  was  put  into  use  recently; 
thus,  it  is  deeper.  The  dashed  lines  in  Figure  23  give 
the  extent  of  the  sumps  as  extrapolated  from  the  ra¬ 
dar  data. 

A  100-MHz  profile  of  the  longitudinal  line  and 
its  migration  are  shown  in  Figure  24.  The  profile  was 
recorded  at  500-ns  time  range  in  order  to  determine 
the  depths  to  the  sewage  in  the  new  sump  and 
to  the  brine  layer  beneath  it.  The  migration  ve¬ 
locity  corresponds  to  an  nm  of  1.46  (e  =  2.13),  but 
nm  values  as  large  as  1.68  also  produced  diffrac¬ 
tions  that  match  some  of  tire  hyperbolas  so  that 
there  is  dielectric  variability.  Continuous  plow¬ 
ing  over  many  years  has  removed  tire  lower 
density  snow.  The  success  of  the  migration  is 
apparent  from  the  linearity  of  the  migrated  fea¬ 
tures  and  the  general  absence  of  diffraction  rem¬ 
nants.  The  conspicuous  features  in  the  migra¬ 
tion  are  the  upper  surfaces  of  the  brine  layer  at 
about  30  m  depth,  the  upper  surface  of  contami¬ 
nated  snow  in  the  newer  (deeper,  less  filled) 
sump  at  about  17  to  19  m  depth,  a  (known) 
wooden  box  at  the  newer  sump  at  about  5.5  m 
depth  (sewer  outfall  head),  and  the  sewage  level 
in  the  older  sump  at  about  8  to  10  m  depth.  The 
sewage  in  the  newer  sump  appears  to  extend  up 
to  45  m  laterally  from  its  discharge  point  (center 
line  symbol),  while  sewage  in  the  older  sump 
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Figure  24.  Sewage  area  longitudinal  100-MHz  profile  (bottom)  and  its  migration  using  e  =  2.8. 


appears  to  extend  laterally  about  30  m.  Taped  dis¬ 
tances  at  each  discharge  point  gave  surface-to-liquid 
depths  of  16.8  and  7.2  m  for  the  new  and  old  sumps 
respectively.  Radar  data  very  near  those  locations 
give  depths  of  17  and  7.8  m.  The  small  discrepancy 
for  the  new  sump  is  easily  within  the  change  of  sur¬ 
face  elevation  near  the  new  shaft  because  of  fre¬ 
quent  surface  regrading  in  that  area  during  the  few 
weeks  between  when  the  distances  were  measured 
and  the  GPR  survey  was  conducted. 

A  hot  point  drill  was  used  to  melt  two  holes 
where  shown  in  Figures  23  and  24.  When  contami¬ 
nated  or  brine-soaked  snow  is  encountered,  the 
rate  of  penetration  of  the  heater  slows  considerably. 
Using  this  technique,  we  determined  the  depth  to 
brine-soaked  snow  at  point  and  sewage-soaked 
snow  at  point  P2  to  be  30  and  18.8  m  respectively. 
The  GPR  data  determined  these  point  depths  to  be 
31  and  20  m,  the  discrepancy  attributable  to  a  less 
dense  overburden  than  nearer  the  shaft. 


Some  of  the  reflecting  horizons  in  Figure  24  ap¬ 
pear  to  block  penetration  to  the  layers  below,  and 
we  assume  this  to  be  the  case  for  the  interpreted 
sewage  area  in  Figure  23.  Most  obvious  is  the  block¬ 
age  of  radio  wave  energy  penetration  to  the  brine 
layer  by  the  sewage.  There  is  also  apparent  blockage 
of  penetration  to  sewage  in  the  newer  sump  by  the 
older  sewage  between  72-  and  90-m  distance.  Along 
this  section,  there  appears  to  be  an  extra  time  delay 
to  the  new  sump  and  this  would  be  ascribable  to  the 
higher  e  of  the  sewage  above.  There  is  also  a  near¬ 
surface  reflecting  layer  from  96  to  150  m,  most  likely 
related  to  a  change  in  snow  density,  that  may  pro¬ 
hibit  penetration  to  the  east  side  of  the  older  sump. 

Figure  25  shows  the  migrated  cross  section  pro¬ 
files,  each  of  which  passed  by  a  sewage  outfall  shaft. 
The  cross  section  of  the  older  sump  shows  it  to  ex¬ 
tend  more  to  the  southeast.  The  cross  section  over 
the  newer  pit  also  shows  extension  to  the  south, 
where  the  level  drops  to  almost  20  m  depth. 


Distance  (m) 
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Figure  25.  Migrated,  100-MHz  cross-sectional  profiles  across  the 
new  and  the  old  sewage  sumps. 


SUMMARY  AND  CONCLUSIONS 

Low  density  layers  at  the  Pegasus  runway,  brine 
discontinuities  along  the  Pegasus  access  road,  me¬ 
tallic  debris  at  the  abandoned  solid  waste  dump, 
and  contaminated  snow  levels  within  sumps  and 
their  extent  at  Williams  Field  were  clearly  shown  by 
reflections  and  diffractions.  The  contrast  between 
the  horizontal  bubble  layering  at  Pegasus  and  the 
folded  density  structure  beneath  parts  of  the  access 
road  is  consistent  with  the  ice  and  snow  contrast 
between  the  two  areas.  No  diffractions  were  seen 
along  the  access  road  that  might  indicate  a  major 
disturbance  such  as  crevasses  within  the  shelf  struc¬ 
ture.  The  radar  results  at  the  fuel  spill  site  could  not 
be  confirmed  because  there  has  been  no  drilling  or 
excavation  to  verify  the  findings. 

Of  particular  interest  is  the  extensive  structure 
of  low  density  layering  seen  beneath  the  Pegasus 
runway,  generally  throughout  the  9-m  depth  to 


the  brine  intrusion.  Previously,  only  lim¬ 
ited  coring  has  revealed  the  presence  of 
these  bubble  layers,  but  the  radar  profiles 
show  how  pervasive  they  are.  Each  of  the 
profiles  in  Figures  8-10  are  at  successively 
greater  resolution  and  reveal  successively 
more  detail,  but  the  time  and  distance 
scales  of  Figure  9  seem  best  suited  for  re¬ 
vealing  the  overall  structure  of  thin  layers 
extending  several  tens  of  meters.  A  deep 
core  and  thin  section  analysis  to  at  least  9 
m  depth  would  be  helpful  in  determining 
the  evolution  of  this  structure  and  the  na¬ 
ture  of  the  ice  or  snow  that  allows  brine  to 
penetrate  below  the  runway. 

Diffractions  from  discontinuities  in  brine 
penetration  within  the  ice  shelf,  as  well  as 
from  the  sewage  deposits,  responded  well 
to  the  single-layer  migration  algorithm 
despite  obvious  variability  in  dielectric 
structure.  The  success  of  the  migrations 
is  judged  by  the  absence  of  coherent  arti¬ 
facts  (concave  or  convex  hyperbolas)  that 
characteristically  occur  when  the  migra¬ 
tion  velocity  is  too  high,  too  low,  or  layer¬ 
ing  with  significant  velocity  contrasts  is 
present.  This  means  that  the  refractive  in¬ 
dex  values  nm  used  were  effective  aver¬ 
ages  for  the  overlying  snow  and  ice. 

The  nm  values  for  the  upper  10  m  at  the 
Williams  Field  sites  range  from  1.38  to 
1.52.  This  range  agrees  with  the  average  n 
of  about  1.4  for  the  top  10  m  or  so  of  the  ice 
shelf  derived  from  the  1978  vertical  den¬ 
sity  profiles  in  this  area  given  by  Kovacs  et  al.  (1982). 
The  nm  values  for  snow  beneath  the  Pegasus  access 
road  are  slightly  higher.  Only  the  upper  meter  or  so 
of  snow  along  this  road  has  been  compacted,  so  that 
a  core  analysis  would  be  needed  to  verify  if  densities 
in  this  area  are  greater  than  along  the  line  investi¬ 
gated  by  Kovacs  et  al. 

The  low  loss  of  the  snow  and  ice  media  should 
allow  bandwidths  centered  as  high  as  1000  MHz  to 
further  improve  vertical  resolution.  Such  transduc¬ 
ers  are  commercially  available  as  are  data  densities 
of  2048  and  4096  samples/scan  and  16-bit  quantiza¬ 
tion  to  allow  such  bandwidths  to  acquire  quality 
data  at  time  ranges  in  excess  of  200  ns.  The  Polar 
Operations  Section,  Division  of  Polar  Programs, 
National  Science  Foundation,  should  consider  pur¬ 
chasing  a  GPR  system  for  use  in  Antarctica  as  these 
studies  in  the  McMurdo  area  and  others  at  the  South 
Pole  (Arcone  et  al.,  in  prep.2)  indicate  that  it  could  be 
used  to  answer  many  practical  questions. 
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We  also  recommend  that  dielectric  properties  of 
mixtures  of  fuel  and  dense  snow  be  measured  as  a 
function  of  temperature  to  determine  if  such  prop¬ 
erties  differ  significantly  from  those  of  dense  snow 
alone.  All  forms  of  oil  and  fuel  have  very  low  dielec¬ 
tric  constants  and  may  not  be  readily  detectable 
where  they  are  dispersed  in  snow  unless  hydrocar¬ 
bon  properties  change  significantly  with  decreasing 
temperature.  If  concentration  rather  than  dispersal 
occurs,  then  fuel  and  oil  detection  is  routine  because 
far  more  subtle  dielectric  contrasts  are  readily  ap¬ 
parent,  as  was  seen  along  the  road  from  Pegasus 
runway  to  Williams  Field. 
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APPENDIX  A:  AIRBORNE  PROFILE  OF  A  PORTION  OF  THE  ACCESS  ROAD 


During  January  1994,  extensive  surveying  was 
done  in  the  McMurdo  area  to  demonstrate  crevasse 
detection  by  airborne  radar  (Delaney  and  Arcone 
1994).  One  of  the  flights  was  almost  directly  over 
section  C-WIL  along  the  access  road  from  the  Pe¬ 
gasus  runway  to  Williams  Field  (Fig.  12)  and  the 
radar  profile  is  shown  in  Figure  Al.  The  center  fre¬ 
quency  of  the  wavelet  is  approximately  200  MHz. 
The  profile  has  been  corrected  for  flight  altitude 
variations,  deconvolved  and  horizontally  filtered  to 
remove  aircraft  reflections. 


The  airborne  profile  shows  the  main  features 
seen  in  Figure  17,  including  the  rising  brine  level 
between  points  C  and  D  and  at  least  one  of  the  dis¬ 
continuities.  Although  the  airborne  profile  is  appar¬ 
ently  closer  to  the  sea  ice  and  does  not  extend  as  far 
northeast  as  does  the  surface  profile,  the  deepening 
brine  level  beyond  the  discontinuity  near  point  WIL 
seems  to  be  deeper  and  extend  further  than  when 
seen  from  the  surface.  This  is  attributable  to  inaccu¬ 
rate  coordinate  control,  movement  of  the  ice  shelf 
during  the  year  and,  probably,  complicated  brine  in¬ 
trusion  and  ice  shelf  dynamics. 
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Figure  Al.  Airborne  radar  profile  at  200  MHz  near  the  Pegasus  runway  access  road  approximately  be¬ 
tween  points  C  and  WIL  (Fig.  12).  The  depth  scale  is  based  on  e  =  2.7.  The  profile  is  to  be  compared  with  that  in 
Figure  17.  The  profile  has  been  corrected  for  aircraft  altitude  fluctuations  and  has  been  deconvolved.  The  weak 
response  of  the  brine  layer  in  the  southwest  half  may  be  attributable  to  a  rough  surface. 
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